
function Cr3bp_ScnEMRot(aux)
%
% 画CR3BP场景，地月旋转系
%
% 时间：2024年06月03日
% 作者：张晨
% 单位：CSU
% 邮箱：chenzhang@csu.ac.cn
%%%%%%%%%%%%%%%%%%%%%%%

% ---------------- 载入参数 ----------------
mu = aux.EMRot.mu;

% ----------------------- fig1 --------------------------------------
figure1;

% 地球位置
xe = - mu;
ye = 0;
text(xe , ye , 'Earth' , 'linewidth' , 1);

% 月球位置
xm = 1 - mu;
ym = 0;
text(xm , ym , 'Moon' , 'linewidth' , 1);

% 画地球
earthReq = 6378 / aux.dim.EMRot_l;
radius = earthReq;
center = [xe, ye]; % 圆心位置
theta = linspace(0, 2*pi, 100); % 构造圆的角度序列
x = center(1) + radius * cos(theta); % x坐标
y = center(2) + radius * sin(theta); % y坐标
fill(x, y, [0 , 0.4470 , 0.7410]);

% 画月球
moonReq = 1738 / aux.dim.EMRot_l;
radius = moonReq;
center = [xm, ym]; % 圆心位置
theta = linspace(0, 2*pi, 100); % 构造圆的角度序列
x = center(1) + radius * cos(theta); % x坐标
y = center(2) + radius * sin(theta); % y坐标
fill(x, y, [0.5 , 0.5 , 0.5]);

% 画LEO和LLO
leoRmag = aux.dim.rMagLEO / aux.dim.EMRot_l;
plot_o([xe , ye] , leoRmag , 'k--' , 0.5);

lloRmag = aux.dim.rMagLLO / aux.dim.EMRot_l;
plot_o([xm , ym] , lloRmag , 'k--' , 0.5);

% 画平动点
[Li_pos , ~] = Cr3bp_Li(mu);
plot(Li_pos(1,1) , Li_pos(1,2) , 'r+' , 'linewidth' , 1.2);
plot(Li_pos(2,1) , Li_pos(2,2) , 'r+' , 'linewidth' , 1.2);
plot(Li_pos(3,1) , Li_pos(3,2) , 'r+' , 'linewidth' , 1.2);
plot(Li_pos(4,1) , Li_pos(4,2) , 'r+' , 'linewidth' , 1.2);
plot(Li_pos(5,1) , Li_pos(5,2) , 'r+' , 'linewidth' , 1.2);

text(Li_pos(1,1) , Li_pos(1,2) - 0.001, 'L1' , 'linewidth' , 1.2);
text(Li_pos(2,1) , Li_pos(2,2) - 0.001, 'L2' , 'linewidth' , 1.2);
text(Li_pos(3,1) , Li_pos(3,2) - 0.001, 'L3' , 'linewidth' , 1.2);
text(Li_pos(4,1) , Li_pos(4,2) - 0.001, 'L4' , 'linewidth' , 1.2);
text(Li_pos(5,1) , Li_pos(5,2) - 0.001, 'L5' , 'linewidth' , 1.2);

xlabel('x / LU');
ylabel('y / LU');
view(0 , 90);
title('Earth-Moon Rotating')

axisBd = 50e4 / aux.dim.EMRot_l;
axis([- axisBd , axisBd , -axisBd , axisBd]);

end
